#/usr/bin/Rscript

on_100 <- read.table("on/profits.dat")
op_100 <- read.table("rf/profits.dat")
ut_100 <- read.table("ut/profits.dat")
iosup_100 <- read.table("iosup/profits.dat")
super_100 <- read.table("super/profits.dat")

print("100 0 on")
shapiro.test(on_100$V1)
print("100 op")
shapiro.test(op_100$V1)
print("100 ut")
shapiro.test(ut_100$V1)
print("100 iosup")
shapiro.test(iosup_100$V1)
print("100 super")
shapiro.test(super_100$V1)

#QQplots

jpeg("normalidade_100.jpg")
par(mfrow=c(2, 3))

qqnorm(on_100$V1, main="on 100 0")
qqline(on_100$V1)
qqnorm(op_100$V1, main="op 100 0")
qqline(op_100$V1)
qqnorm(ut_100$V1, main="ut 100 0")
qqline(ut_100$V1)
qqnorm(iosup_100$V1, main="iosup 100 0")
qqline(iosup_100$V1)
qqnorm(super_100$V1, main="super 100 0")
qqline(super_100$V1)

jpeg("homocedasticidade_100.jpg")
par(mfrow=c(2, 3))

plot(on_100$V1, main="on 100 0")
plot(op_100$V1, main="op 100 0")
plot(ut_100$V1, main="ut 100 0")
plot(iosup_100$V1, main="iosup 100 0")
plot(super_100$V1, main="super 100 0")

dev.off()

